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Abstract 

The detailed analysis of the global structure of the renormalization-group (RG) flow 
diagram for a model with isotropic and cubic interactions is carried out in the frame- 
work of the massive field theory directly in three dimensions (3D) within an assump- 
tion of isotropic exchange. Perturbative expansions for RG functions are calculated 
for arbitrary N up to the four-loop order and resummed by means of the generalized 
Pade-Borel-Leroy technique. Coordinates and stability matrix eigenvalues for the cubic 
fixed point are found under the optimal value of the transformation parameter. Critical 
dimensionality of the model is proved to be equal to N c = 2.89 ± 0.02 that agrees well 
with the estimate obtained on the basis of the five-loop e-expansion [H. Kleinert and V. 
Schulte-Frohlinde, Phys.Lett. B342, 284 (1995)] resummed by the above method. As 
a consequence, the cubic fixed point should be stable in 3D for N > 3, and the critical 
exponents controlling phase transitions in three-dimensional magnets should belong to 
the cubic universality class. The critical behavior of the random Ising model being the 
nontrivial particular case of the cubic model when N = is also investigated. For all 
physical quantities of interest the most accurate numerical estimates with their error 
bounds are obtained. The results achieved in the work are discussed along with the 
predictions given by other theoretical approaches and experimental data. 
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I Introduction 



At present the critical behavior of the basic models of phase transitions described by 
isotropic field theories with a quartic interaction like Y.f=i( L p'i) 2 is we U studied in the 
framework of different theoretical approaches. In particular, the critical phenomena 
in polymers, easy-axis ferromagnets, simple liquids, and binary mixtures, easy-plane 
ferromagnets, certain superconductors, as well as superfluid helium-4, Heisenberg fer- 
romagnets, and quak-gluon plasma in some models of quantum chromodynamics were 
proved to be governed by the O(N) symmetric universality class with N = 0, 1,2,3, 
and 4, respectively. The large-order field-theoretical renormalization group (RG) ex- 
pansions combined with proper resummation techniques, the high-temperature series 
method, and the most advanced Monte Carlo (MC) simulations provided high-precision 
and comparable to each other numerical estimates for important physical quantities 
(critical exponents, universal couplings, and critical amplitude ratios) which nowadays 
are regarded as canonical numbers fl], |2| . 

However in real crystals, due to their complex crystalline structure, an anisotropy 
is always present. That is why, when studying phase transitions in real substances, 
besides the 0(N) symmetric term one should take into account additional quartic 
interactions in corresponding fluctuation Landau- Wilson (LW) Hamiltonian. The sim- 
plest nontrivial crystalline anisotropy is a cubic one. Proper quartic interaction term 
is represented as Y^iLi <pf- In this case the vector of magnetization is directed either 
along the edges or diagonals of a hypercube in N dimensions of the order parameter 
field. 

The critical thermodynamics of magnetic and structural phase transitions in three- 
dimensional (3D) cubic crystals has been extensively investigated more than twenty 
years ago in a good number of papers. By using the lower-order RG approximations, 
Wilson and Fisher ||, Aharony [|J, and Ketley and Wallace || ||] showed that in the 
critical region the fluctuation instability of continuous phase transitions may be ob- 
served, and that it may lead to the isotropization of the system with a cubic anisotropy. 
This fact gave rise to a question what regime of the critical behavior is actually realized 
in 3D cubic crystal with N = 3. Much efforts of many people have been devoted to 
answer this question. It was understood soon that for the given model it is enough 
to calculate the so-called critical (marginal) dimensionality iV c of the order parameter 
field. Indeed, the critical value N c separates two different regimes of critical behavior 
of the system. For N > N c the cubic rather than the isotropic fixed point is stable in 
3D. At iV = N c the points interchange their stability so that for N < N c the stable 
fixed point is the isotropic one. Therefore the calculation of N c is the crucial point in 
studying the critical phenomena in 3D cubic crystal. However attempts to evaluate 
the critical dimensionality resulted in the dramatically different estimates. 
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In fact, the field-theoretical RG analysis of the stability matrix eigenvalues of the 
cubic and isotropic fixed points fulfilled in the one-loop approximation as well as some 
symmetry arguments (for details see Sec. Ill of the paper) lead to the conclusion that 
N c should lie between 2 and 4. Many years ago the three-loop expansion for N c as a 
power series in e has been obtained in Ref . [f| . Summation of that short series at e = 1 
(D = 3) by means of the Pade approximant [1/1] yielded the value N c = 3.128 0, while 
making use of the Pade-Borel resummation method results in the estimate N c = 3. In 
contrast to this, in the work Ref. ||, by using the variational modification of the 
Wilson recursion relation method, it has been found that N c = 2.3. Later, however, 
Newman and Riedel by decoupling the infinite system of the recursion relations for the 
so-called scaling fields and then solving them showed || that for D = 3 N c ~ 3.4. At 
the same time, the classical technique of the high-temperature expansions, under some 
circumstances, allowed to establish that for N = 3 the isotropic critical asymptotics 
in the cubic crystal is unstable [10], thus implying that iV c < 3. Further, ten years 



ago the analysis of the critical behavior of the (mn)-component field model, which has 
a good number of interesting applications to the phase transitions in real substances, 
has been carried out within the three-loop RG approach in three-dimensions. The 
calculation of the stability matrix eigenvalues for the cubic model (m = 1, n = 3) 
provided the stability of the cubic fixed point in 3D, and the critical dimensionality 
turned out to be equal to 2.91 [11]]. In agreement with this, the estimate N c = 2.9 was 
quoted in the work Ref. ||12|| . More recently, Kleinert and Schulte-Frohlinde calculated 
the RG functions for the cubic model in (4 — e) dimensions up to the five-loop order 
H|. Resummation of the critical dimensionality expansion with the help of the Pade 



approximant [2/2] gave the estimate N c = 2.958 [14]. The cubic fixed point eigenvalues 
found by means of a simple resummation algorithm of the Borel type, accounting the 
large-order behavior of the /^-functions when the parameter of anisotropy is very small 



[|iq[ , indicated that the cubic point is stable in 3D fL? |. Finally, in the recent work Ref. 
[ |D| by using finite size scaling techniques and the high-precision MC simulations it 
has been suggested that N c coincides with three exactly. So, such strong scattering in 
the estimates of N c inspired us to study this problem with particular care. Calculation 
of the critical dimensionality of the order parameter field as well as the eigenvalue 
exponents for the cubic fixed point by exploiting the higher-order RG approach in 
three dimensions and generalized Pade-Borel-Leroy (PBL) resummation technique is 
the main goal of the paper. As will be shown, our estimates for N c and eigenvalues are 
in excellent agreement with recent results by Kleinert and collaborators [I7|] obtained 
on the basis of the five-loop e-expansions. 

The paper is organized in the following way. In Sec. II the perturbative expansions 
for /3-functions of the hypercubic model are deduced within the RG technique in 3D 
up to the four-loop order. In Sec. Ill the structure of the RG flows of the model is 
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investigated and fixed point locations for N > 3 are calculated using the generalized 
PBL resummation method. The eigenvalue exponents of the most intriguing O(N)- 
symmetric and cubic fixed points are evaluated for the physically important case N = 3 
and their stability problem is analyzed in Sec. IV. The numerical estimate of the critical 
dimensionality N c , at which the topology of the flow diagram changes, is obtained by 
resumming both the four-loop RG expansions for the /9-functions in 3D and the five- 
loop e-expansion for N c at e = 1. In Sec. V the four- loop RG expansions of the 
critical exponents for generic N and their numerical estimates are presented. Sec. VI 
is devoted to the study of the critical behavior of the three-dimensional random Ising 
model (RIM) which is the special case of the cubic model when N = 0. The coordinates, 
eigenvalues and the critical exponents of the RIM fixed point are computed by applying 
the PBL resummation procedure. The correction-to-scaling exponent uo is estimated 
therein. The results of the investigation are discussed in Conclusion, along with the 
predictions and numerical estimates obtained earlier on the basis of the same or other 
theoretical approaches and experimental data. 



II The model and /3-functions 

The fluctuation LW Hamiltonian of the model reads 

( r 1 1 

H = J d z x[-(ml$ + d^ipid^i) + —(u Gl jkl + u G^)</Wi</Wz 

where <fi, i = 1, . . . , N is the real vector order parameter field in 3D and is the 
linear measure of the temperature, uq and vq denote the "bare" coupling constants. 
The symmetrized tensors associated with isotropic and cubic interactions are 

respectively. When the cubic symmetry is present, the anisotropic exchange has been 



shown to be negligible within the e-expantion method | 19| . We assume that the 
anisotropic exchange is an irrelevant variable in our case too, and study the critical 
behavior of model ([D with isotropic exchange only. 

The model (|l|) has a number of interesting applications to the phase transitions in 
simple and complicated systems. In fact, when N = 1 Hamiltonian ([!]) describes the 
critical phenomena in pure spin systems (the pure Ising model), while for N = 2 it 
corresponds to the anisotropic XY model (the model of two coupled scalar fields) de- 
scribing as structural phase transitions in ferroelectrics as ordering the two-component 
alloys H The magnetic and structural phase transitions in a cubic crystal are 
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governed by model ([!]) as N = 3. Further, when N = Hamiltonian ([!]) determines 
the critical properties of weakly disordered quenched systems undergoing second-order 
phase transitions. The latter is the nontrivial specific case of the hypercubic model, 
the systematical studying of which was initiated in the classical works by Harris and 
Lubensky 



21| [22j and Khmelnitskii |23fl and then considerably advanced by many au- 
thors when employing the conventional field-theoretical RG approach both in 3D and 
(4 — e) dimensions. Finally, the case N — > oo corresponds to the Ising model with 
equilibrium magnetic impurities [24]. In this limit the Ising critical exponent of spe- 
cific heat a changes its sign and takes the Fisher renormalization as well as v and 
7: a — > — Oi/{l — a), v — > v/(l — a), 7 — > 7/(1 — a). Since the critical phenomena in 
the pure Ising and anisotropic XY models are well understood, we will consider below 
the critical behavior of the cubic and random Ising models only. 

To calculate RG functions the standard normalization conditions of the massive 
renormalized theory at fixed dimensions are applied []26|] . For each Feynman graph 
contributing to RG functions the corresponding contractions are computed by using 
the algorithm developed in Ref . |27] . The combinatorial factors as well as the integral 
values are known from Ref. p8[. After performing simple but cumbersome calculations 



we obtain the four-loop expansions for the /3-functions: 



0n 



U < 1 — U 



6 



N 



v + 



1 



3 (2.024691 iV + 9.382716) 



3 (0.449648 



u 



+ 44.444444 u v + 10.222222 v' 

(iV + 8) 3 

+ 18.313459 N + 66.546806) u 3 + 3 (6.646878 N + 164.613849) u 2 v 
+ 3 (0.621889 N + 100.955929) u v 2 + 65.937285 v 3 

+ — 1 [-(0.155646 N 3 - 35.820204 N 2 - 602.521231 N 
[N + 8) 4 L 

- 1832.206732) u A - 3 (1.352882 N 2 - 182.073890 N 

- 2064.170701) u 3 v + 3 (27.250336 N + 2110.408809) u 2 v 2 

+ 9 (1.291017 N + 308.599361) u v 3 + 495.005747 v 4 



(3) 



Pv 



V 1 - 



1 



A^ + 8 



'12 u + 9 v) + 



(N + 



+ 92.444444 u v + 34.222222 



1 



(N + 



(3.407407 N + 54.814815) u 2 
-(1.251107 A^ 2 



41.853902 N - 469.333970) u 3 + 9 (0.248784 N + 136.511768) u 2 v 



4 



+ 957.781662 u v 2 + 255.929737 v 3 + — - g ^ 4 (0.574653 

- 0.267107 N 2 + 584.287672 N + 5032.692260) u 4 + 3 (0.057375 N 2 
+ 107.641680 N + 5989.283536) u 3 v + 3 (7321.464604 

- 16.494003 N) u 2 v 2 + 11856.956858 u v 3 + 2470.392521 v 4 



(4) 



It is reasonable to note that cubic model ([!]) possesses some specific symmetry 
property when N = 2. Namely, the transformation of the field components 



(<Pi + <&) 



combined with substitution of the quartic coupling constants 



u — > u + - V , V 



(5) 



(6) 



does not change the structure of the initial Hamiltonian itself. As a result, /5-functions 
(|3[) and (|j) should obey certain symmetry relations [2Ti] which may be written down as 



/ 3 \ 3 
(3 u (u + -v,-vj = (3 u (u, v) + - f3 v (u,v) 

[3 v (u+-v,-v^ = -(3 v (u,v) . 



(7) 



It can be easily verified that equations (0) are really satisfied. The specific symmetry 
of Hamiltonian (HD will be used below to obtain the lower boundary value of the critical 
dimensionality N c . 



Ill Method of summation and fixed point locations 

Before analyzing the four-loop /3-functions let us consider briefly the predictions fol- 
lowing from the one-loop approximation. It is easy to see that /5-functions @ and @ 
in this approximation have four different solutions corresponding to Gaussian (trivial), 
Ising, isotropic (Heisenberg) and cubic fixed points with the coordinates: 
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respectively. The most intriguing fixed points are the isotropic and cubic ones. At 
N = 4 the coordinates of these two points coincide that leads to the conclusion that 
the critical dimensionality N c has the upper boundary value N c = 4. On the other 
hand, the eigenvalue exponents for the isotropic and cubic fixed points are given by 
the expressions: 

Af = -1, A? = (A-4)/(A + 8) , 

Af = -1 , A^ = (4 - iV) / 3N . (9) 

If the real parts of both eigenvalues are negative, the corresponding fixed point is 
infrared stable; if eigenvalues are of opposite sings, the point is of "saddle-knot" type. 
It is seen from Eqs.([|) that N c = 4 actually separates two different regimes of critical 
behavior of the model. When N > N c the cubic rather than the isotropic fixed point 
is stable in 3D, while for N < N c the stable fixed point is the isotropic one. 

To determine the lower boundary of N c one should employ the above mentioned 
specific symmetry property of the model, when N = 2 []30 |. As was already pointed 



out, the rotation @ of the components of <pi by 7r/4 combined with substitution (f|) 
generates relations (0), but does not change the form of the RG equations. However 
for N = 2 transformations (|J), (^) result in the relocation of the coupling constants 
values so that the cubic and Ising fixed points are transformed into each another at the 
3D RG flow diagram. Since the exact RG equations always have the Ising fixed point, 
which inevitably is the saddle-knot one, these equations should have also the cubic 
fixed point, which will be unstable. In this situation, the isotropic fixed point, again 
always existing in the exact RG equations, should be the stable knot only. Therefore 
we conclude that the lower boundary of iV c is not less than two. Of course, the real 
value of N c can be obtained only on the basis of the thorough analysis of the structure 
of the RG flow diagram, provided that the /3-functions of the model are calculated in 
sufficiently high-order RG approximations and then processed by means of appropriate 
resummation techniques. 

Let us now concentrate our attention on the analysis of the four-loop /^-functions. 
It is well known that field-theoretical RG expansions are divergent. The character of 
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their large-order asymptotic behavior for the case of simple 0(iV)-symmetric models 



was established in Refs. [[IT], [32], |33] . In particular, it was proved that the coefficients of 
the series at large k behave as c{—a) k\k , where the asymptotic parameters a, b and c 
are assumed to be calculated for each RG function. Knowledge of the exact values of the 
asymptotic parameters in combination with the most powerful resummation procedure 



of the Borel transformation with a conformal mapping ||34j| , first proposed in Ref. |35] 



and then elaborated in Refs. [gg, |37], [38fl , made it possible to develop the accomplished 



quantitative theory of critical behavior of simple systems [1], 0, |35|, |38|, |39], [40], [41] 



At the same time, the asymptotic nature of RG functions of anisotropic models is 
unknown. Calculation of the exact values of the asymptotic parameters characterizing 
the large-order behavior of the series in such models is a very difficult and still unsolved 
problem. As an exception one should mention the anisotropic quartic quantum oscil- 
lator representing one- dimensional ip A field theory with a cubic anisotropy. Within an 
assumption of the weak anisotropy, the transformation parameters for the perturbative 
expansion of the ground state energy of this system as well as for the /3-functions of 



the cubic model have recently been found 0, |42| . Later this information has been 
used to solve the stability problem of the cubic fixed point in three dimensions from 
e-expansion |17| . 



Usually, in a lack of any information about the high-order behavior of the series 
either the simple Pade-Borel or Chisholm-Borel resummation procedures are used, for 
treating the perturbative expansions of anisotropic models. The latter technique, how- 
ever, possesses at least two inherent drawbacks. First, some ambiguity in the calcula- 
tion of coefficients of denominators of the Chisholm approximants is unavoidable [[43] . 
Second, the Chisholm-Borel procedure does not hold the specific symmetry properties 
of a model. At the same time, exploiting the Borel transformation in combination 
with the Pade or Chisholm approximants shows that the results of calculation are very 
sensitive to the choice of the type of approximants. This may lead to the estimates, 
which do not provide reliable predictions even in the higher-loop RG approximations 
(see, for instance, Sec. VI of the paper). Besides, in the framework of both schemes it 
is very difficult to determine any error bounds for evaluated quantities. 

In the present work we apply for processing RG expansions of the /3-functions 
and critical exponents of model ([I]) the PBL resummation method generalized for the 
two coupling constant case. This resummation technique, first introduced by Baker, 
Nickel, and Meiron in Ref. H , turned out to be highly efficient when used to study 
the critical behavior of simple 0(iV)-symmetric models in 3D. The critical exponent 
estimates obtained within the framework of this technique are regarded nowadays as 
the most accurate values, as those of Ref. |2|, We motivate our choice of the 

PBL resummation method by the following reasons: 

• 3D RG expansions for the /^-functions and critical exponents of the cubic model 
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are alternating in signs. Therefore using the PBL resummation technique is quite 
natural. 



• It can be expected that for complex models with more than one coupling constant, 
the asymptotics of RG series at large orders will comprise a factor k\k b . The PBL 
resummation method removes divergences of this type. 

• The PBL resummation method allows one to determine the error bounds for the 
physical quantities to be calculated, in a natural way. 

The generalized PBL resummation procedure consists of the following steps. Let a 
physical quantity F(u, v) be represented by a double series 

F(u,v)=J2fou> i v j , (io) 

where coefficients ~ {i + + j) b at large orders (i, j — > oo), the additional 
parameter b being an arbitrary non-negative number to be defined below. Associated 
with the initial series ( |T0D is the function 



T{u,v;b) = J e-H b B(ut,vt)dt (11) 
o 



The Borel-Leroy transform B(x,y) is the analytical continuation of its Taylor series 

B(x, y) = £ — - /; ; , /V (12) 
Y T(« + j +b+ 1) 

absolutely convergent in a circle of the nonzero radius. In order to calculate the integral 
in (|TTD one should continue analytically B(x,y) for < x < oo and < y < oo. To 
this end, the rational Pade approximants [L/M] (x, y) are used. The Pade approximant 
method is determined in a conventional way ||15|| . Let us consider a "resolvent" series 

B(x, y, A) = V \ k V /If"' / : = V A ^ k , (13) 
1 ^ &T(A; + 6 + l) to 

where coefficients Ak are uniform polynomials of kth order in u and v. The sum of the 
series is then approximated by 

B(x,y) = [L/M] I (14) 
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The Pade approximants [L/M] in A are given by an attitude 

lL/M] ~ S ' (15) 

where -Pl(A) and Qm(A) are polynomials of degrees L and M respectively with coeffi- 
cients depending on x and y which should be determined from the conditions: 

Q M (X)B(x, y; A) - P L (A) = 0(A L+M+1 ) , 

Qm (0) = 1 . (16) 

Replacing variables x = ut and y — vt in the Pade approximants and then evaluating 
the Borel-Leroy integral 

oo 

F(u, v;b)= J e~H b [L/M] | x _dt . (17) 
o 

we obtain the approximate expressions for RG functions. 

Among the Pade approximants the diagonal (L = M) or near-diagonal ones were 
proved to exhibit the best approximating properties (T^|. However, as the degree of 
the denominator M increases, the number of possible poles of the approximant in- 
creases too. If some of the poles belong to the positive real semiaxis, the corresponding 
approximant should be rejected. Due to this, the choice of "working" approximants, 
which might be used for analytical continuation of the Borel-Leroy image onto the 
complex cut plane, is largely limited. On the other hand, varying the free parameter b 
in Borel-Leroy transformation ([0]) allows one to optimize the resummation procedure 
under the condition that the fastest convergence of the iteration process is achieved. 
So, taking into account the above mentioned remarks, in order to find the fixed points 
locations of the model we adopt the following scheme. For the fixed N, the /3-functions 



are resummed by virtue of transformation ( IT ) in the highest-loop orders by shift- 
ing the transformation parameter b. In order to make an analytical continuation of 
the Borel-Leroy transforms B u (u,v), B v (u,v) over the cut-plain, the most appropri- 
ate Pade approximants [2/1], [3/1], and [2/2] are chosen. The fixed points locations 
are then determined for each b from solution of a set of equations: /^ es (w c , t> c ) = 0, 
j3l es {u c , v c ) = 0. The final locations are obtained by averaging over the values given by 
the approximants under the optimal value of the parameter b, at which the quantity 
1 1 — J-"l(u, V] 6)/JF l _ 1 (m, v; b) \ reaches its local minima. The quantity Tl{u, v\ b) is eval- 
uated for the L-partial sum of the series in Eq . (|i~7f) , L stands for the step of truncation 
of the series. 

In Fig. 1 the results of computation of the cubic fixed point locations depending 
on the parameter b are presented for the physical important case iV = 3. Three curves 
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correspond to the three Pade approximants. The parameter b shifts from to 3. As 
seen from the figure the optimal value of b is zero. At this point the numerical values 
of the cubic fixed point locations given by different approximants are the most close to 
each other. The result of computing the cubic fixed point locations are also presented 
in Table I. In the first three columns of the table the fixed point locations values found 
for the Pade approximants [2/1], [3/1], and [2/2] at b = are placed. Averaging 
the results of processing over the all approximants under the optimal value of b gives 
the estimates standed in the fourth column of the table. These numbers we adopt 
as the final estimates of the cubic fixed point locations found within the four-loop 
approximation. As an accuracy for these approximate values we take the maximum 
deviations of the average values of the fixed point locations from those given by the 
approximants at b = 0. 

One can observe, looking at Fig. 1, that the values of the cubic fixed point locations 
given by the symmetric approximant [2/2] weakly depend on the shift parameter b. 
Averaging over the all values given by this approximant within the interval [0,3] results 
in the cubic fixed point locations estimates presented in the fifth column of Table I. 
The coordinates of the cubic fixed point found earlier on the basis of the three- and 
four-loop approximations with the use of the Chisholm-Borel resummation method are 
presented in the sixth and seventh columns of the table for comparison. These numbers 
include the normalizing multiplier y needed to compare our /9-functions with those 
obtained in Refs. jnj [Lj|. 

In order to verify the correctness of the chosen approach let us apply the above 
considered scheme to estimate the fixed point locations of the 0(iV)-symmetric model 
where the numerical results are well known. Consider, for example, an 0(3)-symmetric 
case relevant to the Heisenberg ferromagnets. The six-loop 3D RG expansion for the 
/5-function of this model was reported in Refs. [^, The PBL resummation of 
that series with the use of eight types of the Pade approximants [2/1], [3/1], [2/2], 
[4/1], [3/2], [5/1], [4/2], and [3/3] yields, after solving the equation (3 res (g c ) = 0, the 
picture displayed in Fig. 2. It is seen that the values of the isotropic fixed point 
location calculated in the highest RG orders with the help of the approximants [3/3], 
[4/2], and [3/2] are very weakly dependent on the parameter b varied within the interval 
< b < 15. The curves corresponding to these approximants are intersected at b = 4.5. 
Therefore the value b = 4.5 is the optimal value of the transformation parameter in 
which the fastest convergence of the iteration procedure is ensured. For b = 4.5 the 
central value estimate of the isotropic fixed point is g c = 1.392. The maximum deviation 
of the central value from the values given by some of the approximants [3/3], [4/2], 
and [3/2] at the point b = 10 is adopted approximately as an apparent accuracy of the 
calculation, A = 0.0013. Such a small error can be explained by the small dispersion 
of the curves within the range 5 < b < 10. So, the estimate g c = 1.3920 ± 0.0013 is in 
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excellent agreement with those found more then twenty years ago in Refs. [|35], |3| as 
well as with recent results of Ref. [§]. 

Within the framework of the four-loop approximation there are only three appro- 
priate Pade approximants. Averaging the results of computing the isotropic fixed point 
location given by the approximants [2/1], [3/1], and [2/2] under the optimal value of 
the transformation parameter results in the estimate g c = 1.3925 ± 0.0070. The error 
was determined again through the maximum deviation of the central value from those 
given by each of the approximants at b = 0. It is seen that the four-loop estimate of 
the coordinate of the isotropic fixed point is in a good accordance with the best ones 
followed from the six-loop consideration. 

Let us note that the coordinate of the 0(3)-symmetric fixed point calculated within 
the five-loop approximation does not approach the exact value. Namely, the PBL 
resummation procedure leads to the estimate g c = 1.3947 ±0.0040. Although the error 
of the calculation became visibly smaller, the central value of the fixed point location 
stepped aside from the four- and six-loop ones. 

Thus, the fulfilled numerical analysis shows that the isotropic fixed point location 
estimate obtained in the four- loop level occurs to be close to the "exact" value. There- 
fore it can be expected that in the case of the cubic model the fixed point locations 
u c = 1.3428 ± 0.0200, v c = 0.0815 ± 0.0300 (see fourth column of Table I) will not 
distinguish strongly from "exact" values as well. The coordinates of the cubic fixed 
point for some values N of the order parameter dimensionality are presented in Table 
II. Our calculations show that for iV = 3 the coordinates of the cubic fixed point prac- 
tically do not differ from those of the Heisenberg one. However, with increasing N the 
cubic fixed point runs away from the isotropic point moving towards the Ising one. In 
the large N limit these two fixed points become close to each another so much that 
the influence of the 0(iV)-symmetric invariant on the critical thermodynamics of the 
cubic model vanishes. This can be easily seen by applying the consideration to the 
one-loop solutions of the RG equations of model ([[]). Indeed, rescaling the coupling 
constants u — > u/N, v —>■ v/N m the initial Hamiltonian and taking then the limit 
N — ► oo in Eqs.(Bp one can see that the cubic fixed point approaches the Ising one 
asymptotically. So, the cubic model turns out to be splitted into iV non-interacting 
Ising models, the critical behavior of each of them is determined by a set of critical 



exponents renormalized according to Fisher [25 



The data listed in Table II will be used further for calculating the stability matrix 
eigenvalues as well as the critical exponents of the hypercubic fixed point. 
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IV Stability and critical dimensionality 



One of the independent way to determine fixed point locations in fixed D is to construct 
the RG flow phase diagram of a model. If, at the flows diagram, there exists a fixed 
point of stable knot type, the trajectories originated from some point within the range 
of stability of the initial Hamiltonian would flow towards the knot. The region at the 
flow diagram where the trajectories are intersected provides the coordinates of this 
stable fixed point. Thoroughly investigating the 3D RG flow diagram of cubic model 
(|IJ) in the four-loop approximation we arrived at the conclusion that the cubic rather 
than the isotropic fixed point is absolutely stable for all N > 3. 

On the other hand, the reliable conclusion about the stability of the cubic fixed 
point for N > 3 can be given on the basis of calculating the eigenvalue exponents A's 
of the stability matrix 



taken at u = u c and v = v c . If the real parts of both eigenvalues are negative, the fixed 
point is the stable knot on the (u,v) plane. If Ai, A2 have opposite signs, the point is 
of the "saddle-knot" type. 

To calculate the stability matrix eigenvalues of the cubic and isotropic fixed points 
we have chosen the following strategy. First, the derivatives of the /3-functions (0), @ 
are calculated, and new RG expansions resummed by means of the PBL technique are 
substituted into the matrix My. The eigenvalue exponents of the matrix of derivatives 
Mij obtained in this way are evaluated then under the optimal value of the transfor- 
mation parameter b. In Fig. 3 a), b) we present our numerical results for the series 
— and — for the physically interesting case N = 3. The curves correspond to the 
three types of the Pade approximants used within the four-loop approximation. The 
crossing of the curves gives the optimal value of b at which we find ^ \ opt = —0.7536 
and ^f|opt = —0.0331. Because the series — and — are turned out to be shorter 
by one order in comparison with — and — their resumming performed with the 
help of the approximant [2/1] only yields the monotonic dependence of the result of 
processing on the parameter b. In this unfavorable situation, we take into account an 
additional Pade approximant [1/1] to optimize the iteration procedure. The results 
are plotted in Fig. 3 c), d). For the optimal values of b we obtain ^§^\ pt = —0.4566 
and ^fflopt = —0.0409. Straightforward calculation of the eigenvalues of the stability 
matrix gives for the cubic fixed point the numbers placed in Table III. The eigen- 
values of the isotropic fixed point as well as the analogous numerical estimates obtained 
recently in Ref. |I7j on the basis of using the five-loop e-expansions are presented for 
comparison therein. These estimates show that the cubic fixed point is absolutely sta- 
ble in 3D for N = 3, while the isotropic fixed point appears to be stable on the it-axis 
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only. Our numerical results agree well with those obtained in Ref. |17| . 

Let us now calculate the critical dimensionality iV c of the order parameter field. 
The critical dimensionality is defined as a value of N at which the cubic fixed point 
coincides with the isotropic one. Equivalently, for N = N c the second eigenvalue of the 
stability matrix My vanishes, A2 = 0. 

Studying carefully the 3D RG flow diagram of model ([|) depending on the order of 
approximation with the use of different Pade approximants we arrive at the conclusion 
that N c = 2.910 ± 0.035 and N c = 2.890 ± 0.020 within the three- and four-loop 
approximations respectively. The accuracy of calculation of N c was determined through 
the evaluation of the stability matrix eigenvalues for different N from the interval of the 
above mentioned errors. That value of N = N c , above or below of its central number, 
at which the second eigenvalue A2 was becoming nonzero, was taking for the upper or 
lower boundary of N c respectively. 

It is worthy to compare the four-loop estimate of N c just found with that which can 
be obtained within the e-expansion method. The five-loop e-expansion for N c has been 
calculated in Ref. ||13[| . The series turned out to be alternating in signs that allows one 



to resum it by means of the PBL technique. To this end, we will use again the most 
appropriate Pade approximants [2/1], [3/1] and [2/2] for analytical continuation of the 
Borel-Leroy transform for all < et < 00. Dependence of the results of processing 
the critical dimensionality N c on the transformation parameter b is depicted in Fig. 4. 
The curves corresponding to the approximants are crossed at the point b ~ 1. The 
appropriate value of the critical dimensionality is N c = 2.894 ± 0.040. As an error 
of the calculation it is natural to assume the maximum scattering of numerical values 
given by the approximants at b = with respect to the value obtained at the crossing 
point of the curves. This estimate of N c is in excellent agreement with the above 
found within the 3D RG approach. So, both schemes, the RG technique directly in 3D 
and the ^-expansion method, result in the same estimate of the critical dimensionality 
N c = 2.89, thus implying that the cubic fixed point is stable in three dimensions for 
N > 3. This means that the critical behavior of model (|l|) should be governed by the 
cubic fixed point with a certain set of critical exponents which will be calculated for 
generic N in the next section. 



V Critical exponents for generic N 

Having the coordinates of the cubic fixed point, the stability of which in 3D for N > 3 
has been proved in the previous section, it is possible to obtain the numerical estimates 
for the critical exponents. Two of them are known to determine the critical behavior 
of the system while the others can be found via the famous scaling laws. The four- 
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loop RG expressions for the magnetic susceptibility exponent 7 1 as well as for the 
correlation function exponent 77 for generic symmetry index N are as follows: 



7 



-1 



1 



iV + 8 
1 



(N + 2) 3 

u + - v 

2 2 
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0.024684 (N + 2) (N + 8) u 3 + 0.222156 (N + 8) u 2 v 

1 



+ 1.99940 u v z + 0.666468 v e 



+ 



(N + 8Y 



(19) 



(0.004299 A^ - 0.667985 A^ 



- 4.60922 N - 6.51210) u 4 - 3 (0.017194 A^ 2 - 2.70633 N - 13.0242 ) u 3 v 
+ 3 (0.681615 N + 22.8884) u 2 v 2 + 47.140 u v 3 + 11.7850 v A 



These series, however, are known to be divergent. To extract from them a physical 
information concerning the critical behavior of the substances of interest we will apply 
the same resummation procedure as that used above. Note, since the coefficients of 
the series of r] are rapidly diminishing, its values are found by the direct substitution 
of the coordinates of the cubic fixed point into Eq. ([l9|) , whereas the series of the 
exponent / y~ 1 before the substitution needs to be resummed. Using the generalized 
PBL resummation technique, we evaluate the magnetic susceptibility exponent 7 and 
then, having the values of 77 at hand, estimate the correlation length critical exponent 
v by the scaling relation. The result of numerical processing Eq. ( |I8|) depending on the 
parameter b for A^ = 3 is dipicted in Fig. 5. The three types of Pade approximants 
have been used, and the transformation parameter shifted within the range [0,15]. For 
the optimal value of b we obtain 7 = 1.3775 ± 0.0040. Taking into account that in the 
four-loop approximation rj = 0.0332 ± 0.0030, we arrive at the estimate of the critical 
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exponent v: v — 0.6996±0.0037. The error bounds for the exponents 7 and v have been 
determined through a maximum deviation of their central values at b = b opt from the 
corresponding values given by the approximants at b = 0. In the case of the exponent 
T], as the uncertainty we take the absolute difference between two successive, three- and 
four-loop, results. The critical exponents estimates of the cubic fixed point for other 



values of N are listed in Table II. The estimates obtained from the e-expansions |45 
are placed for comparison therein. 

How much do the critical exponent estimates found differ from the "exact" ones? To 
answer this we resort again to the analysis of numerical results for the 0(3)-symmetric 
model obtained on the basis of using the PBL procedure. In Fig. 6 the results of 
processing the 3D RG series for the susceptibility exponent 7 in successive orders of 
perturbation theory in a number of loops are presented. Within the four-loop approx- 
imation under the optimal value of transformation parameter we find 7 = 1.3778. The 
six-loop calculations under the optimal choice of b (the point of crossing the curves 
corresponding to the approximants [4/2], [5/1], and [4/1]) give 7 = 1.3867. The latter 
estimate seems to be in a good agreement with that known from Refs. [[35|, lifl . Thus, 
the difference between the four-loop results and their six-loop counterparts does not 
exceed 0.01. 

For N = 2 the numerical results for the cubic model turn out to be much better. 
In this case the cubic fixed point lies in the v < half plane and is unstable (N < N c ). 
However for N = 2 the Hamiltonian of model (|I|) possesses the specific symmetry [see 
Eqs.(^|), (H)] that transforms the cubic point into the Ising one and vice versa. Due 
to the symmetry, the critical exponents of both fixed points coincide. In the four- 
loop level for the cubic fixed point at iV = 2 we found the estimates: 7 = 1.2416, 
7] = 0.0323, and v = 0.631. These numbers agree well with the best estimates followed 



from the six-loop RG expansions P, |35|, |3^, [0]]. One may hope therefore that the 



critical exponent estimates for the cubic model obtained in the present work within 
the four- loop approximation will differ from the "exact" values by no more than 1-2%. 

Unfortunately, for the most important case N = 3 having a good number of inter- 
esting applications to the critical phenomena in real substances, the critical exponents 
of the cubic and isotropic fixed points turn out to be practically the same. This is 
the consequence of the closeness of both points on the 3D RG flow diagram. Thus, 
although in the course of the investigation the cubic fixed point has been shown to be 
stable at iV = 3 and, therefore, the critical behavior of the magnetic phase transitions 
in crystals with cubic anisotropy should belong to the cubic rather than the isotropic 
universality class, a certain difficulty arises in trying to identify the cubic fixed point 
from experimentally determined exponents. Due to this "near-marginality" , the cal- 
culation of the critical exponents in cubic magnets seems to be of academic interest 
only. 
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At the same time, the numerical analysis shows that as N increases the distinction 
between the critical exponents of the cubic and isotropic fixed points increases as well. 
In the limit iV — > oo the critical exponents of the cubic fixed point go over into those 
of the Ising model with equilibrium magnetic impurities |25| . 



VI Numerical results for the RIM 

In the critical region the character of critical behavior of the defect crystals such as 
impure uniaxial ferromagnets LiTbi- x Y x F4, Cdi^ x Nd x Cl^ or diluted Ising antiferro- 
magnets Mn x Zrii_ x F2, Fe x Zn\- x F2 are known to be described by the 3D random 
Ising model with the effective Hamiltonian 

H = J rfx[i(m^ 2 + (V^) 2 ) + ^ 4 + <A(xV 2 ] , (20) 

where ^(x) is the static random field describing fluctuations of local transition tem- 
perature, m\ — m\ c ~ T — T c . Averaging over all configurations of with Gaussian 



weight and employing the replica trick |46| one can reduce problem ( f20|) to the analysis 
of critical behavior of the iV-component hypercubic model ([I]) in the limit N —>■ [[J7| . 
Moreover, the Ising vertex vq in Eq.(^) plays a role of the cubic vertex in Eq.([I]), 
while a role of the impure vertex plays the isotropic one uq. Obviously, because now 
Wo < 0,t>o > 0, the added interaction of critical fluctuations of the order parameter 
field through the intermediary of impurities is the attraction. 

Studying the magnetic and structural phase transitions in weakly disordered sys- 
tems is of considerable interest both from theoretical and experimental point of view. 
It is well known that the critical exponents of such systems should differ markedly from 



those of the pure ones, due to the famous Harris criterion |48[ . Indeed, according to 
Harris, provided the specific-heat exponent a purc of a pure system is positive, i.e. the 
specific-heat is divergent at the critical point (C ~ \T — T c \~ a ), a new critical behavior 
under dilution is expected. Otherwise, if a purc < 0, the critical behavior of the ran- 
dom system would be similar to that of the pure one. Among the three-dimensional 
iV-vector models, only the Ising model has a > 0, and the corresponding new critical 
behavior has been obtained as an unusual stable fixed point dependent upon y/e, where 
e = 4 — D. Furthermore, one (smallest) of the eigenvalues of the stability matrix taken 
at the RIM fixed point is expressed through the critical exponents a and v: A2 ~ oc/u 
(for the details see H |5Q|). 

The systematical calculation of the RIM critical exponents was hystorically begun 
in the seminal works by Harris and Lubensky |2l], 22 and Khmelnitskii |[23|| . However 



the e-expansion technique could not provide the reliable numerical estimates, because 
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RG equations of model ([I]) for N = turn out to be degenerate in the one-loop 
approximation. Such a degeneracy causes powers of y/e to appear in expansions for the 
fixed point locations as well as for the critical exponents, thus leading to the substantial 
decrease of accuracy expected within the high-loop approximations |51], [52[ . 

The following pronounced step to evaluate the RIM exponents was made by G.Jug 



}I| , who applied the alternative approach, the RG in fixed dimensions. The reasonable 



numerical estimates were obtained within the two-loop approximation by making use of 
the Chisholm-Borel procedure to resum the resulting series. Later the critical exponents 
series for the 3D RIM were deduced within the three- and four-loop approximations 
and corresponding numerical estimates were obtained on the basis of the Chisholm- 



Borel summation method |TT], Pade-Borel procedure and the first confluent form 
of the e algorithm of Wynn |M| . 

In this section we will study the critical thermodynamics of the 3D RIM using 
the generalized PBL resummation method. Setting N = in Eqs. (^) and (^) and 
solving the system of equations /3£ es (w c , v c ) = 0, (3l es {u c , v c ) = we find the RIM fixed 
point locations depending on the transformation parameter b. Resulting curves are 
depicted in Fig. 7. Unfortunately, this picture is not complete, because in the three- 
loop approximation the Pade approximants [2/1] and [1/2] have the poles for all b. On 
the other hand, the fixed point locations values given by the approximant [1/1] seem to 
be very far from the true ones. So, in order to determine the RIM fixed point locations 
we cannot apply the optimization algorithm described in Sec. Ill, at least within the 
given approximation. In such situation we need a new working criterion. We can, for 
instance, select the approximants providing the most stable values under the variation 
of the parameter b. As is seen from Fig. 7 the locations of the RIM fixed point given 
by the approximants [2/2], for the w c -component, and [3/1], for the t> c -component, are 
practically independent on the parameter b. Indeed, the dispersion of the corresponding 
curves within the range < b < 4 is no more than 1 ■ 10 -4 . The fixed point locations 
obtained in such a way may then be used for calculation of the critical exponents 7, 
rj, and v as well as of the stability matrix eigenvalues when starting the optimization 
procedure (see Sees. IV and V). Corresponding numerical estimates are summarized in 
Table IV. For comparison we collected in the table the data obtained earlier either by 
resumming RG functions within the minimal subtraction scheme directly at D = 3 (3D 
MS) f55| , or by applying different resummation procedures (Chisholm-Borel technique 
[fill , [H, 0, £ algorithm of Wynn, "AW" @), or found experimentally ]57], |58| and 



through the MC simulations [|59|]. To compare those results with our numbers, the 
fixed point locations are given by taking the normalizing factor | into consideration. 
As experimental data we consider the averaged values of Ref. |57| obtained in a course 
of studying the critical behavior of the site-random Ising system Mn x Zn\_ x F2 with 
x = 0.75 (or 0.50) by the neutron scattering method as well as the averaged values of 
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Ref. |5|. 

In Table IV we present also the numerical results for the fixed point locations and 
critical exponents obtained on the basis of the simple Borel summation method com- 
bined with the Pade approximants [3/1] and [2/2]. Although the fixed point coordinates 
obtained in such a way turned out to be strongly different, the critical exponents esti- 
mates differ nevertheless from each other only slightly. From this point of view neither 
approximant is better. 

The another possible way to determine the coordinates of the RIM fixed point is to 
calculate them as averages between the values given by the highest approximants [3/1] 
and [2/2] for each of the components (see Fig. 7). The corresponding critical exponents 
estimates are found to agree with just considered as well as with experimental data 
and MC results, within the error bounds (see Table IV). 

Unlike the cubic fixed point, the determination of the error bounds for the RIM 
fixed point locations is a more difficult problem. In the case of selecting the most 
stable approximants we have taken the following scheme. First, the values of the fixed 
point locations given by the approximants [3/1] and [2/2] are averaged for each of 
the components separately within the interval < b < 4. The discrepancy between 
averages in these approximants is then adopted as a sought uncertainty in the results. 
At the same time, the error bounds for the critical exponents are determined in the 
same way as in Sec. V. If we adopt as the fixed point locations the averages between 
the values given by the highest approximants [3/1] and [2/2] separately for each of the 
components, the error bounds seem to be even smaller, they are almost three quarters 
of the previous ones. Absence of the error bounds in some of the places of Table IV 
means that the errors either cannot be established or were not established. 

We have checked also the stability of the RIM fixed point on the three-dimensional 
RG flow diagram. In all considered cases the eigenvalues of the stability matrix turned 
out to be negative, except the calculations based on the simple Borel summation 
method with the Pade approximant [2/2]. In this case the second eigenvalue A2 oc- 
curred to be positive and too large. This is in contradiction to the known theoretical 
and experimental predictions as well as to the Monte-Carlo simulations Indeed, 
the second, smallest in modulos, eigenvalue A2 of the matrix of derivatives of the f3- 
functions is well known to define the so called correction-to-scaling exponent u that 
governs the leading corrections to the universal power laws. Thus, the approach of the 
zero field susceptibility to the critical temperature, for T > T c , is characterized by the 
Wegner series |60| 

x~r T-f(i + r 1 7*»' + ...), (21) 

with Tfc being the nonuniversal amplitudes and r = (T — T c )/T c |61| . Here 7, v are the 
asymptotic values of the susceptibility and correlation length critical exponents. As the 
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exponent u decreases the region increases, where the corrections to scaling laws should 
be taken into account. So, the smallness of uj in the RIM indicates the importance of 
its calculation for analysis of the asymptotic critical behavior of dilute systems . 

Recent MC calculations based on the analysis of the first correction term in Eq. (|2ip 
provided the estimate of the correction-to-scaling exponent uj = 0.37±0.06 |59| . Almost 
the same number was more recently obtained in the framework of the four-loop 3D RG 
analysis used for processing divergent series the Borel summation method in combina- 
tion with the simple rational Chisholm approximants like [M,M/1,1], uj = 0.372±0.005 
[|6^| . Although the apparent accuracy of this estimate seems to be highly overstated, 
the central value is in accordance with previous estimates uj = 0.366 p| and uj = 0.359 



11JI, u = 0.376 |56[ derived within the three-loop approximation in the framework of 
the minimal subtraction scheme and the 3D RG respectively. Our estimate of uj ob- 
tained on the basis of the Borel summation method with the Pade approximant [3/1] 
is close to the above mentioned. On the contrary, using the Pade approximant [2/2] 
in the Borel transformation leads to the unphysical result for the correction-to-scaling 
exponent. 

At the same time, applying the PBL resummation method to study the asymptotic 
critical behavior of the systems with impurities results in the correction-to-scaling 
exponent values which are different from those predicted by either MC simulations or 
simple resummation procedures (see Table IV). It is the consequence of the fact that 
the four-loop approximation is not enough to obtain reliable estimates of the RIM fixed 
point locations. Note, however, that our estimates of uj are within the error interval 
found for the MC result. The critical exponents estimates obtained confirm also the 



inequality — vuj < a < conjectured for the random models |B3|. Unfortunately, at 
present we cannot indicate any error bounds in our calculation of the exponent uj. 

So, if to assume that the MC simulations |59| provide the numerical estimate of uj 
which is close to the "exact" one, a question is to be put forward: can the estimation 
of the correction-to-scaling exponent of the dilute systems be used as an additional 
criterion of selection of the resummation techniques to be employed? Probably the 
answer will be given in the course of the further investigation of the critical properties 
of the RIM within the higher-order RG approximations provided a more sophisticated 
method of the series summation will be used, on a level with the simple techniques. 
Thus more recently, a new approach to summation of divergent field-theoretical series 
has been suggested p5 |. The method, based on the Borel transformation combined 



with a conformal mapping, relies upon the stability of the result of processing on the 
transformation parameters and therefore does not require knowing the exact asymptotic 
behavior of the series. This method has been tested on the functions expanded in 
their asymptotic power series and applied to estimating the ground state energy of 
quantum mechanical systems, including anisotropic oscillator, as well as to calculating 
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the critical exponents for some conformal field theories |45|, |65| . The successful testing 
the developed technique on simple systems made it possible to apply it for obtaining 
the critical exponents estimates of anisotropic iV-vector field models describing both 
magnetic and structural phase transitions in cubic and tetragonal crystals |27|, f45f . It 



would be reasonable therefore to apply this technique in studying the critical behavior 
of the RIM exactly in three dimensions within the five-or higher-loop approximations. 



Conclusion 

The complete RG analysis of a field model with two quartic coupling constants associ- 
ated with isotropic and cubic interactions describing magnetic and structural phase 
transitions in a good number of real substances has been carried out within the 
four-loop approximation directly in three dimensions. Perturbative expansions for 
/3-functions and critical exponents were deduced for generic N. The fixed points loca- 
tions were found for N > 3 by applying the generalized Pade-Borel-Leroy resummation 
technique, and the global structure of the 3D RG flow diagram was investigated. The 
analysis of the eigenvalue exponents of the most intriguing isotropic and cubic fixed 
points fulfilled for the physically important case N = 3 has shown that the cubic rather 
than the isotropic fixed point is absolutely stable in 3D. The eigenvalues estimates of 
both fixed points were found to agree well with those of recently calculated on the basis 



of exploiting the five-loop e-expansions |T7[. The critical dimensionality N c of the order 



parameter field, at which the topology of the flow diagram changes, has been analyzed 
by the two different methods: 1) by resumming the four- loop RG expansions for the 
/3-functions in 3D and 2) by resumming the five-loop e-expansion for N c at e = 1. 
The numerical estimates N c = 2.89 ± 0.02 and N c = 2.894 ± 0.040 obtained are in a 
good agreement with the earlier results |Tl|, [12|] and confirm the conclusion about the 
stability of the cubic fixed point for iV > 3. Consequently, the magnetic and structural 
phase transitions in three-dimensional anisotropic crystals with cubic symmetry are of 
second order and their critical thermodynamics should be governed by the cubic fixed 
point with a specific set of critical exponents. The corresponding four-loop critical 
exponents estimates were found in the framework of the PBL resummation method. 
On the basis of comparative numerical analysis with the 0(iV)-symmetric model, the 
critical exponents of which are solid established |]], it has been shown that in the case 
of the cubic model the difference between the four-loop estimates and the exact val- 
ues does not exceed 1-2%. Although our results for the most interesting case N = 3 
are in good accordance with earlier theoretical predictions, the cubic universality class 
is not easy to distinguish experimentally from the isotropic one, due to the obvious 
marginality of the problem, N c ~ 3. 
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The critical behavior of weakly quenched disordered systems undergoing second 
order phase transitions and described by the three-dimensional random Ising model, 
which is the nontrivial specific case of the cubic model when N = 0, has also been 
investigated. The coordinates, eigenvalues and critical exponents of the RIM fixed point 
were computed by using the PBL resummation method and the more simple Pade-Borel 
procedure. Our numerical results along with the known theoretical and experimental 
data were summarized in the table. Although the RIM fixed point locations found in 
the framework of different approximation schemes turned out to be strongly different, 
the critical exponents estimates differ from each other only slightly, within the error 
bounds obtained. 

A special attention was given to study the stability of the 3D RIM fixed point on the 
RG flow diagram. The calculation of the stability matrix eigenvalues based on applying 
different resummation techniques showed that they are negative. Consequently, the 
RIM fixed pont is stable in 3D. As an exception we have indicated the case of using 
the simple Borel summation method in combination with the Pade approximant [2/2], 
which led to the unphysical result for the second eigenvalue A2 = 00. Applying the 
PBL procedure to calculate the leading correction-to-scaling exponent uj for the three- 
dimensional impure systems was shown to result in the numerical estimates which 
are distinguishable markedly from those predicted by recent MC simulations [59] or 



followed from applying simple Borel-like resummation procedures [XT], [56], |62[ (see Table 
IV). Note, however, that our estimates of to were found to be within the error bounds 
known for the MC result. 

At last, it is worth noting that one can hardly hope to extract a reliable numer- 
ical estimate of u from the five-loop e-expansion even resummed by a proper 
method. The point is that, the e-series for u turns out to be very short, due to the 
degeneracy of the random Ising model /5-functions on the one-loop level. Therefore, 
until the appreciable discrepancy between the results given by different methods of the 
series summation on the one hand, and the MC calculations, on the other, does exist, 
the further investigation of the asymptotic critical behavior of dilute systems will be 
highly desirable in the framework of the higher-order (five-or six-loop) RG approxi- 
mations provided a more sophisticated resummation procedure, for instance the Borel 
transformation combined with a conformal mapping, will be applied. 

Note added to manuscript. More recently the six-loop study of critical behavior of 



the 3D cubic model appeared |66| . Perturbative expansions for /^-functions and critical 
exponents deduced within the massive field theory in fixed dimension have been re- 
sumed by means of the Borel transformation combined with a conformal mapping that 
takes into account the singularities of the Borel transform. The fixed point locations, 
stability matrix eigenvalues and critical exponents estimates obtained turned out to 
be essentially the same results as those of the present work. In Appendix we present 
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our numerical estimates for the cubic model obtained on the basis of the six-loop ex- 
pansions of Ref. 
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when the PBL resummation procedure is applied. The critical 
behavior of three-dimensional random Ising systems has recently been studied within 
the five- and six-loop RG approximations in Ref. |67| and Ref. p8| respectively. 



Acknowledgments 

It is pleasure to thank Drs. A. I. Mudrov, B. N. Shalaev and Prof. A. I. Sokolov for 
a reading of the manuscript. I am especially grateful to Dr. B. N. Shalaev for many 
illuminating discussions on properties of the phase transitions of disordered systems. I 
thank Dr. E. Vicari for useful correspondence and comments. 



Appendix 

In this section we present our numerical estimates for the cubic fixed point locations 
and critical exponents for some values of N (see Table V) using the six-loop expansions 
for RG functions recently obtained by J. M. Carmona, A. Pelissetto, and E. Vicari 
pfifl . We apply the generalized PBL resummation technique. The susceptibility and 
correlation length critical exponents are estimated through the original series for 7 -1 
and v, whereas the critical exponent 77 is obtained by the known scaling relation: 
rj = 2 — 7/V. It is seen that our estimates are in excellent agreement with those of Ref. 

where the Borel resummation procedure with a conformal mapping, that takes 
into account the singularities of the Borel transform, has been applied. 



References 

[1] J. Zinn- Justin Quantum Field Theory and Critical Phenomena, Clarendon Press 
(Oxford 1989, third ed. 1996) and related articles in this book. 

[2] R. Guida and J. Zinn- Justin, J. Phys. A31, 8103 (1998). 

[3] K. G. Wilson and M. E. Fisher, Phys. Rev. Lett. 28, 240 (1972); K. G. Wilson, 
Phys. Rev. Lett. 28, 548 (1972); 

[4] A. Aharony, Phys. Rev. B8, 4270 (1973). 

[5] I. J. Ketley and D. J. Wallace, J. Phys. A6, 1667 (1973). 



22 



[6] D. J. Wallace, J. Phys. C6, 1390 (1973). See also R. A. Cowley and A. D. Bruce, 
J. Phys. C6, L191 (1973). 

[7] A. Aharony, A. D. Bruce, Phys. Rev. Lett. 33, 427 (1974); A. Aharony (unpub- 
lished) . 

[8] M. C. Yalabik and A. Houghton, Phys. Lett. A61, 1 (1977). 

[9] K. E. Newman and E. K. Riedel, Phys. Rev. B25, 264 (1982). 
[10] M. Ferer, J. P. Van Dyke, and W. J. Camp, Phys. Rev. B23, 2367 (1981). 
[11] N. A. Shpot, Phys. Lett. A142, 474 (1989). 

[12] I. O. Mayer, A. I. Sokolov, and B. N. Shalaev, Ferroelectrics. 95, 93 (1989). 

In this article the four-loop RG expansions of the model are reported for the cases 
N = and N = 3 only. The critical exponents estimates for the cubic fixed point 
have been obtained in three dimensions by making use of the Chisholm-Borel 
summation method. However somewhat analysis for N c as well as for eigenvalue 
exponents of the cubic fixed point has not been given. 



[13 
[14 

[15 

[16 
[17 
[18 
[19 
[20 
[21 
[22 
[23 



H. Kleinert and V. Schulte-Frohlinde, Phys. Lett. B342, 284 (1995). 

It was established that the diagonal or near-diagonal Pade approximants exhibit 
the best approximating properties. See, for instance, |y| 



G. A. Baker, Jn. and P. Graves-Morris, Pade Approximants (Reading, MA: 
Addison-Wesley) 1981. 

H. Kleinert and S. Thorns, Phys. Rev. D52, 5926 (1995). 

H. Kleinert, S. Thorns, and V. Schulte-Frohlinde, Phys. Rev. B56, 14428 (1997). 
M. Caselle and M. Hasenbusch, J. Phys. A31, 4603 (1998). 
A. D. Bruce, J. Phys. C7, 2089 (1974). 

E. Brezin, J. C. Le Guillou, and J. Zinn- Justin, Phys. Rev. BIO, 892 (1974). 
A. B. Harris and T. C. Lubensky, Phys. Rev. Lett. 33, 1540 (1974). 
T. C. Lubensky, Phys. Rev. Bll, 3573 (1975). 

D. E. Khmelnitskii, Zh. Eksp. Teor. Fiz. 68, 1960 (1975) [Sov. Phys. JETP 41, 
981 (1976)]. 



23 



[24] A. Aharony, Phys. Rev. Lett. 31, 1494 (1973). 
[25] M. E. Fisher, Phys. Rev. 176, 257 (1968). 

[26] G. Parisi, in: Proceedings of the Cargrese Summer School 1973 (unpublished); J. 
Stat. Phys. 23, 49 (1980). 

[27] A. I. Mudrov and K. B. Varnashev, Phys. Rev. B57, 3562 (1998). 

[28] B. G. Nickel, D. I. Meiron, and G. A. Baker, Jr., University of Guelph report, 
1977 (unpublished). 

[29] A. L. Korzhenevskii, Zh. Eksp. Teor. Fiz. 71, 1434 (1976) [Sov. Phys. JETP 44, 
751 (1976)]. 

[30] I. O. Mayer, A. I. Sokolov, Izv. Akad. Nauk SSSR, Ser. Fiz. 51, 2103 (1987). 

[31] L. N. Lipatov, Zh. Eksp. Teor. Fiz. 72, 411 (1977) [Sov. Phys. JETP 45, 216 
(1977)]. 

[32] E. Brezin, J. C. Le Guillou, and J. Zinn- Justin, Phys. Rev. D15, 1544 (1977). 
[33] E. Brezin, G. Parisi, J. Stat. Phys. 19, 269 (1978). 

[34] J. J. Loeffel, Centre d' Etudes Nucleaires de Saclay Report No. SACLAY-DPh- 
T/76-20 (unpublished). 

[35] J. C. Le Guillou and J. Zinn- Justin, Phys. Rev. Lett. 39, 95 (1977). 

[36] R. Seznec and J. Zinn- Justin, J. Math. Phys. 20, 1398 (1979). 

[37] D. I. Kazakov, O. V. Tarasov, D. V. Shirkov, Theor. Math. Phys. 38, 15 (1979). 

[38] A. A. Vladimirov, D. I. Kazakov, O. V. Tarasov, Zh. Eksp. Teor. Fiz. 77, 1035 
(1979) [Sov. Phys. JETP 50, 521 (1979)]. 

[39] J. C. Le Guillou and J. Zinn- Justin, Phys. Rev. B21, 3976 (1980). 

[40] S. G. Gorishny, S. A. Larin, and F. V. Tkachov, Phys. Lett. A101, 120 (1984). In 
the article six of the 135 diagrams contributing to the /3-functions and anomalous 
dimensions of the 0(A)-symmetric y9 4 -theory in (4 — e) dimensions were calculated 
erroneously. Their analytic recalculation was performed in Ref. H. Kleinert, J. 
Neu, V. Schulte-Frohlinde, K. G. Chetyrkin, and S. A. Larin, Phys. Lett. B272, 
39 (1991); ibid. B319, 545(E) (1993). However the critical exponent estimates 
found after reevaluation turned out to be very small different from the previous 
ones. 



24 



[41] J. C. Le Guillou and J. Zinn- Justin, J. Rhys. Lett. (Paris) 46, L137 (1985); J. 
Phys. (Paris) 48, 19 (1987); ibid. 50, 1365 (1989). 

[42] H. Kleinert, S. Thorns and W. Janke, Phys. Rev. A55, 915 (1997). 

[43] J. S. R. Chisholm, Math. Comput. 27, 841 (1973). 

[44] G. A. Baker, Jr., B. G. Nickel, and D. I. Meiron, Phys. Rev. B17, 1365 (1978). 
In the Reference section of the article there is a remark that the extension of the 
Borel method of summation for the first time was studied by E. Leroy in Ann. 
Fac. Sci. Univer. Toulouse 2, 317 (1900). 

[45] A. I. Mudrov and K. B. Varnashev, Phys. Rev. E58, 5371 (1998). 

[46] K. Ziegler, Nucl. Phys. B285 [FS19], 606 (1987); Nucl. Phys. B344, 499 (1990). 

[47] G. Grinstein, A. Luther, Phys. Rev. B13, 1329 (1976); A. Aharony, Phys. Rev. 
B13, 2092 (1976). 

[48] A. B. Harris, J. Phys. C7, 1671 (1974). 

[49] J. Sak, Phys. Rev. BIO, 3957 (1974). 

[50] S.-k. Ma, Modern Theory of Critical Phenomena (W. A. Benjamin, Inc., London, 
1976). 

[51] B. N. Shalaev, Zh. Eksp. Teor. Fiz. 73, 2301 (1977) [Sov. Phys. JETP 46, 1204 
(1977)]; C. Jayaprakash, H. J. Katz, Phys. Rev. B16, 3987 (1977); B. N. Shalaev, 
S. A. Antonenko, and A. I. Sokolov, Phys. Lett. A230, 105 (1997). 

[52] A similar phenomenon has recently been observed in the course of studying another 
complicated model having three independent coupling constants and describing 
magnetic and structural phase transitions in a good number of antiferromagnets 
and ferroelectrics. See A. I. Mudrov and K. B. Varnashev, Phys. Rev. B57, 5704 
(1998). Half- integer powers in e arising in that model have a different origin but 
also result in the loss of accuracy. 

[53] G. Jug, Phys. Rev. B27, 609 (1983). 

[54] I. O. Mayer, J. Phys. A22: Math. Gen, 2815 (1989). 

[55] R. Folk, Yu. Holovatch, T. Yavors'kii, J. Phys. Stud. 2, 213 (1998). 

[56] Yu. Holovatch, T. Yavors'kii, J. Stat. Phys. 92, 785 (1998). 



25 



[57] P. W. Mitchell et al, Phys. Rev. B34, 4719 (1986). 



[58] T. R. Thurston, C. J. Peters, R. J. Birgeneau, Phys. Rev. B37, 9559 (1988). 
H. G. Ballesteros et al, Phys. Rev. B58, 2740 (1998). 



[60] F. J. Wegner, Phys. Rev. B5, 4529 (1972); ibid. B6, 1891 (1972) and Critical 
Phenomena and the renormalization group, Lectures of the VIII Finnish Summer 
School (1973). 



[61] Strictly speaking, Eq.(^Tj) was obtained by F. J. Wegner for the pure system. It 
can be shown however, by solving formally RG equations for the cubic model 
when N = 0, that Eq.(^) will be true for the RIM likewise with the correction- 
to-scaling exponent of the pure model replacing the impure one u = A2. For 
the details see E. Brezin, J. C. Le Guillou, and J. Zinn- Justin. Field theoretical 
approach to critical phenomena, in: Phase Transitions and Critical Phenomena, 
edited by C. Domb and M. S. Green. Vol. 6 (Academic Press, New York 1976). 

[62] R. Folk, Yu. Holovatch, T. Yavors'kii, Pis'ma Zh. Eksp. Toer. Fiz. 69, 698 (1999) 
[JETP Lett. 69, 747 (1999)]. 

[63] H. K. Janssen, K Oerding, and E. Sengespeick, J. Phys. A28, 6073 (1995). 

[64] A. Aharony, A. B. Harris, and S. Wiseman, Phys. Rev. Lett. 81, 252 (1998) 

A. I. Mudrov and K. B. Varnashev. "New approach to summation of divergent 
field-theoretical series in models with strong coupling" , in: Problems of Quantum 
Field Theory. (JINR Publishing Department, E2-99-35, Dubna, Russia. Eds.: B. 
M. Barbashov, G. V. Efimov, and A. V. Efremov), 267 (1999); |hep-th/981lT25 . 



[66] J. M. Carmona, A. Pelissetto, and E. Vicari. "The TV-component Ginzburg-Landau 
Hamiltonian with cubic anisotropy: a six-loop study", Phys. Rev. B (2000) in 
press; |cond-mat / 99 1 2 1 L5 . 

[67] D. V. Pakhnin and A. I. Sokolov. "Five-loop renormalizat ion-group expansions 
for the three-dimensional n-vector cubic model and critical exponents for impure 
Ising systems" , |cond-mat/9912071[ 

[68] A. Pelissetto and E. Vicari. "Randomly dilute spin models: a six-loop field- 
theoretic study", |cond-mat/ 0002404 



26 



[69] The stability matrix eigenvalues estimates obtained in the framework of the PBL 
resummation technique with the use of approximants [5/1] and [4/1] for the series 
, and approximants [4/1] and [3/1] for the series , under the optimal 
value of the transformation parameter b are Ai = —0.7887 and A2 = —0.0126 
for N = 3, Ai = -0.7895 and A 2 = -0.0760 for N = 4, A x = -0.8004 and 
A 2 = -0.1414 for iV = 8, Ai = -0.8054 and A 2 = -0.1769 for N = 00. 



27 



Figure Captions 



Fig. 1. Curves demonstrating dependence of the results of calculating the cubic fixed 
point locations a) -u c -component and b) t> c -component on transformation parameter b 
from the four-loop approximation for N = 3. The upper curve (O) corresponds to the 
[2/1] Pade approximant (three-loop approximation), while the middle (A) and lower 
(□) curves correspond to the [2/2] and [3/1] Pade approximants (four-loop approxima- 
tion) respectively. 

Fig. 2. The results of computation of the 0(3)-symmetric fixed point locations from the 
three- to the six-loop approximations obtained on the basis of the PBL resummation 
method with eight types of the approximants: [2/1] - O, [3/1] - □, [2/2] - A, [4/1] - 
full O, [3/2] - full A, [5/1] - full o, [4/2] - o, [3/3] - x. On the six-loop level under 
the optimal value of the transformation parameter b = 4.5 we obtain the estimate 
9c = 1-392. 

Fig. 3. Graphs of dependence of the results of processing the series a) — fff , b) — , 
c) — ^f, and d) — on the parameter b for iV = 3. The curves on the pictures 
a) and b) are given in the same notations as in the previous figures, while for the 
curves corresponding to the approximants [1/1] and [2/1] on the pictures c) and d) the 
notations O and □ are used respectively. 

Fig. 4. Dependence of the results of processing the e-series for the critical dimension- 
ality N c on transformation parameter b. The crossing of the curves corresponding to 
the three Pade approximants gives the optimal value of the parameter b at which the 
estimate N c = 2.894 is obtained. 

Fig. 5. Curves demonstrating dependence of the result of processing susceptibility 
exponent 7 of the cubic model on the parameter b for N = 3 in the four-loop approxi- 
mation. 

Fig. 6. Graphs of dependence of the results of processing the 3D RG series for sus- 
ceptibility exponent 7 of the 0(3)-symmetric model on the parameter b for the eight 
types of Pade approximants in the framework of the PBL procedure. 

Fig. 7. Graphs of dependence of the results of calculating the RIM fixed point locations 
a) -Uc-component and b) t> c -component on transformation parameter b within the four- 
loop approximation. The curves marked by O, A, and □ correspond to the Pade 
approximants [1/1], [2/2], and [3/1] respectively. 
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Table I: Coordinates of the cubic fixed point of RG equations for N = 3 found under 
the optimal value of the transformation parameter 6 = 0. "AV" denotes the average 
value. 





[2/1] 


[3/1] 


[2/2] 


AV 


AV over [2/2] 


Ref. [11] 


Ref. [12| 


u c 


1.3536 


1.3338 


1.3410 


1.3428 
± 0.0200 


1.3425 


1.348 


1.3357 




0.0526 


0.1026 


0.0894 


0.0815 
± 0.0300 


0.0937 


0.090 


0.0906 



Table II: Coordinates of the cubic fixed point of RG equations and critical exponents 
estimates for some iV found under the optimal value of the transformation parameter b 
within the four-loop approximation. The critical exponents values obtained within the 
framework of the e-expansion method (five-loop results) and marked by the symbol (*) 
are presented for comparison. 



N 


u c 


Vc 


V 


V 


7 


3 


1.3428 


0.0815 


0.0332 


0.6996 


1.3775 








0.0375* 


0.6997* 


1.3746* 


4 


0.9055 


0.8167 


0.0327 


0.7131 


1.4028 








0.0365* 


0.7225* 


1.4208* 


5 


0.6980 


1.2361 


0.0325 


0.7154 


1.4076 








0.0358* 


0.7290* 


1.4305* 


6 


0.5807 


1.5386 


0.0324 


0.7157 


1.4082 








0.0354* 


0.7301* 


1.4322* 


7 


0.5060 


1.7874 


0.0324 


0.7155 


1.4079 


8 


0.4544 


2.0076 


0.0324 


0.7153 


1.4074 


9 


0.4168 


2.2108 


0.0324 


0.7149 


1.4067 


10 


0.3881 


2.4032 


0.0324 


0.7147 


1.4063 


12 


0.3473 


2.7680 


0.0323 


0.7142 


1.4054 



* Quoted from Ref. @ 
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Table III: Four-loop eigenvalue exponents estimates for the cubic (CFP) and Heisenberg 
(HFP) fixed points (N = 3) found in 3D under the optimal value of the transformation 
parameter b. 





CFP 


CFP, Ref. ]T7| 


HFP 


HFP, Ref. |T7| 


Ai 


-0.7786 


-0.7648 


-0.7791 


-0.7640 


A 2 


-0.0081 


-0.0085 


0.0077 


0.0089 
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Table IV: Numerical results for the RIM. Here "MSA" indicates that the fixed point 
locations and the critical exponents are found for the most stable approximants in 
the framework of the PBL resummation technique, "AV" means that the fixed point 
locations are calculated as the averages between the values given by the highest approx- 
imants [3/1] and [2/2], "PB" denotes a simple Pade-Borel resummation. The results 
obtained either within the different theoretical approaches, or on the basis of using dif- 
ferent resummation techniques, or experimentally ("Exp.") and by means of the MC 
simulations, are presented for comparison. 





u c 


V c 


V 


V 


7 


OJ 


MSA 


-0.5816 
± 0.0850 


1.9822 
± 0.0740 


0.040 
± 0.011 


0.681 
± 0.012 


1.336 
± 0.020 


0.310 


AV 


-0.6246 
± 0.0600 


1.9438 
± 0.0500 


0.034 
± 0.010 


0.672 
± 0.004 


1.323 
± 0.010 


0.330 


PB [3/1] 


-0.6839 


1.9877 


0.033 


0.674 


1.326 


0.362 


PB [2/2] 


-0.5800 


1.8934 


0.034 


0.669 


1.316 




AW 


-0.5874 a 


1.9362 a 




0.668 a 


1.318 a 




Ref. [11] 


-0.728 


2.006 


0.021 


0.671 


1.328 


0.359 


Ref. @ 


-0.745 


2.011 


0.019 


0.671 


1.328 


0.376 


Ref. [Q 


-0.6668 


1.9951 


0.034 


0.670 


1.326 




3D MS 






0.053" 


0.677" 


1.319" 


0.330" 
0.390 e 


Exp. 








0.71 ± 0.02 c 
0.70 ± 0.02 a ' d 


1.37 ± 0.04 c 
1.37 ± 0.04 a ' d 




MC, 
Ref. [0 








0.6837 
± 0.0053 


1.342 
± 0.010 


0.37 
± 0.06 



a Quoted from Ref. §4] 
b Quoted from Ref. || 
c Quoted from Ref. (57) 
d Quoted from Ref. §§ 
e Quoted from Ref. || 
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Table V: The cubic fixed point locations and critical exponents estimates for some N 
obtained within the PBL resummation procedure in the six-loop approximation. Here 
the coupling constants are given in the same notations as those of Ref. |66 . 



N 


u c 


V c 


V 


V 


7 


Ai [69] 


a 2 m 


3 


1.3177 
± 0.0170 


0.0964 
± 0.0165 


0.0327 
± 0.0020 


0.7040 
± 0.0040 


1.3850 
± 0.0050 


- 0.7833 
± 0.0054 


- 0.0109 
± 0.0032 


4 


0.8804 
± 0.0080 


0.6360 
± 0.0050 


0.0316 
± 0.0025 


0.7150 
± 0.0050 


1.4074 
± 0.0030 


- 0.7887 
± 0.0090 


- 0.0740 
± 0.0065 


8 


0.4410 
± 0.0070 


1.1331 
± 0.0160 


0.0305 
± 0.0025 


0.7143 
± 0.0035 


1.4068 
± 0.0030 


- 0.7955 
± 0.0150 


- 0.1396 
± 0.0100 


oo 


0.1751 
± 0.0040 


1.4122 
± 0.0090 


0.0319 
± 0.0035 


0.7094 
± 0.0030 


1.3962 
± 0.0040 


- 0.7986 
± 0.0200 


- 0.1787 
± 0.0050 
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